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Abstract. A relatively fast algorithm for evaluating Weil-Petersson 
volumes of moduli spaces of complex algebraic curves is proposed. 
On the basis of numerical data, a conjectural large genus asymp- 
totics of the Weil-Petersson volumes is computed. Asymptotic for- 
mulas for the intersection numbers involving f/;-classes are conjec- 
tured as well. The precision of the formulas is high enough to 
believe that they are exact. 



The aim of this note is to report on the recent progress in com- 
puting Weil-Petersson volumes of moduli spaces of complex algebraic 
curves (with or without marked points) that resulted from better pro- 
gramming, software and hardware as compared to [7]. The numerical 
evidence led us to a plausible guess about their large genus asymp- 
totic behavior that may have further applications in algebraic geome- 
try, combinatorics, dynamical systems and string theory. 

I. Algorithms 

Let M. 9t n denote the moduli space of stable n-pointed genus g com- 
plex algebraic curves. The universal curve p : M. g>n +i — > M. g>n has n 
canonical sections Xi,...,x n given by the marked points. Put ipi = 
Ci(x*u) G H 2 (M. gn) Q), where ui is the relative dualizing sheaf on 
JA g>n +i. The first Mumford class of M. g<n is the direct image class 

«1 = P*^n + l = / £ H 2 (M g , n , Q). 

J fiber 

The Weil-Petersson metric is Kahler on M. g>n - Its symplectic form 
ujwp extends to M. g , n as a closed current and represents the class 
2tt 2 k 1 G H 2 (M 9i n, R) (see [S]). By definition, the (normalized) Weil- 
Petersson volume of M. g , n is just its standard symplectic volume with 
respect to the form -^ujwp'- 

' ' „3g-3+n 



For all g,n > with 2g + n > 3 these are positive rational numbers. 
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Below we describe an algorithm for computing Weil-Petersson vol- 
umes (see [7] for details). 



Theorem 1. Let 

t V<9y y dt ) 1 " L dt 
where 

{-If- 1 y k 



x(y) = -y/yfoVy/y) = J2 

k= 



■f (Jfe-l)! k\ 

(Jo denotes the Bessel function of the first kind). 
Then 

(i) the KdV equation 



x 2 12 

has a unique solution u(y, t) = y + h 2g Y^Lx u ghJi where each u g (y, t) 
is a Laurent polynomial in t of the form 

59-1 

u g(y,t) = Yl U 9-Ay)t~ k ; 

k=2g+l 

(ii) for each g > 2 the equation 

dl<t> g (y,t) = u g (y,t) 
has a unique solution of the form 

5<?-5 

<f> 9 {y,t) = Yl ^k{y)t~ k \ 

k=2g-l 

(Hi) for any g,n > the Weil-Petersson volume of M. g n is given by 
the formula 

Vg, n = dQ(f) g (y,t) \ y= o :t= i . 

The above theorem extends to the intersection numbers involving 
■0-classes. Fix a set d = (c?i, . . . , d n ) of non- negative integers and put 
\d\ — d 1 + ■ — h d n . Consider 

T/ _ / j.di , d n 3g-3+n-\d\ (r) s 

V >"« - (3,-3 + n-|i|)l L ' ' ' ' (2) 



Theorem 2. Let 

do = Xi [ ^- + x\ (x 2 + — ) + Y^Xfc+i-^ 

\oy v y / dxi t^> dxk 



di = x\ + yd 

OX i 



LARGE GENUS ASYMPTOTICS OF WEIL-PETERSSON VOLUMES 



with x(y) = — v /yJo(2y / ^) as above. 
Then 

(i) the KdV equation 
has a unique solution 

oo 

v(y, x h x 2 ,...)=y + h 2g ^ v g(y> x u x 2, ■ ■ ■ ), 

9=1 

where each v g (y,t) is a polynomial in x\ of the form 

59-1 

v g (y,x 1} x 2 ,...)= ^2 v g;k(y,x 2 ,x 3 , . . . )ac*; 

fe=2g+l 

(wj /or eac/i g > 2 the equation 

dfygiv, x i, x 2 ,...) = v g (y, x 1: x 2 , . . . ) 
/ias a unique solution of the form 

59-5 

il> g (y,x 1 ,x 2 , ■ ■ ■) = Vtofc(y>z 2 ,a: 3 ,...)zi; 
fe=2flr — a 

(raj t/ie intersection number V g ^ n -d is given by the formula 

V g ,n;d = , , d °-ip g (y, xt, x 2 ,...) \ y=0tXl=ltX2=X3= ... =0 , 

0X1 ox 2 ■ ■ ■ 

where Ik is the number ofdi's equal to k. 

The proof follows the same lines as that of Theorem 1 and utilizes 
an observation of M. Kazarian on how to explicitly express mixed in- 
tersection numbers of ip- and K-classes in terms of intersection numbers 
of ^-classes alone [3J. The details will appear elsewhere. Note that for 
d\ — ■ ■ ■ = d n — it reduces to Theorem 1 with the obvious change of 
variable x\ = 1/t. 

The main advantage of our algorithm is its speed, and in this respect 
it is superior to the algorithms of C. Faber pQ and M. Kazarian [3J, 
though it loses to both of them in generality. 

2. ASYMPTOTICS 

It may be instructive to begin with the large n asymptotics of Weil- 
Petersson volumes. The following exact asymptotic formula was proven 
in [1] for any fixed g: 

V g , n = n\ C n n^- 1 ^ 2 {a g + O (1/n)) , n -> 00, (3) 

where C = —zqJq(zo) and z is the first positive zero of the Bessel 
function Jq(z). The coefficients a g can also be explicitly computed [4 J 
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(in fact, one can even get the complete asymptotic expansion of V g<n as 
n — > oo). 

The problem seems more challenging when n is fixed and g — > oo. 
We implemented the algorithm of Theorem 1 in a Maple™ Q program 
and computed all numbers V 9}n for g < 50 and 1 < n < 4. These data 
led us to 

Conjecture 1. For any fixed n > 



This formula agrees with the earlier results of [21 [5]. Approximate 
values of the constants c n with n < 4 are: Co ~ 1.8, c\ ~ 0.75, C2 ~ 
0.1, c 3 « -0.15, c 4 « -0.001. 

Our Maple™ implementation of the algorithm of Theorem 2 eval- 
uates the intersection numbers V gtn;c t given by ([2]). In particular, we 
computed all V gtn; d with g < 40 and 4 < 2, fc = 1, 2, 3), and from that 
we get 

Conjecture 2. For any fixed n > and a /fzed set d = (di, . . . d n ) o/ 
non-negative integers 



where Ik denotes the number of d t 's that are equal to k. 

Both these conjectures hold numerically with a high precision, so 
there is a good reason to believe that they are actually true. However, 
at the moment these asymptotic formulas lack theoretical justification. 
The next section contains some data and heuristics. 



It is known that the order of magnitude of V^ n is (2g) \ for large g (see 
[21 [5]), so the problem was to find the asymptotics up to the factors 
of smaller order. A question of M. Mirzakhani about the behavior 
of the ratio Vg_i in+ 2/V^ n as g — > oo served us as a starting point. 
Computations show that this ratio decreases with g for n — 0, 1, and 
increases for any fixed n > 2. At the same time for any fixed g this ratio 
decreases when n grows. Below is the table of decimal approximations 
(rounded up to 10 digits) for n = 1,2 and 41 < g < 50: 



©Maplesoft, a division of Waterloo Maple Inc. 
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3. Numerics 
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Vg-1,3/ Vg,l 


T/ , , /!/ „ 
^g-1,4/ ^,2 


A "1 

41 


19.78811999 


19.68915456 


A O 

42 


19.78695294 


19.69037678 


A O 


19.78584026 


19.69154073 


A A 

44 


i A 70 /I 7700 -1 

19.78477824 


ly.oyzo5U4y 


A C 

45 


19.78376350 


in rno Tn n ^ /1 

19.69370974 


4b 


1 n 7007000/1 

ly. / o2 (9294 


in o /i 70 1 
19.694/2186 


47 


19.78186376 


19.69568993 


48 


19.78097335 


19.69661676 


49 


19.78011934 


19.69750494 


50 


19.77929954 


19.69835682 



The behavior of V r 9 _i in+2 /V^ jn suggests that there is a limit as g — > 
oo independent of n. D. Zagier numerically identified this limit with 
2ir 2 = 19.7392088 . . . (private communication). 

The next step is to analyze the behavior of the ratio 2gV^ jn _i/V^ n : 
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2gv g ,i/v g>2 


2gV 9)2 /V 9t3 


2gV g>3 /V gA 


40 


0.5082406948 


0.5031382404 


0.4981366818 


41 


0.5080365079 


0.5030613837 


0.4981822417 


42 


0.5078421948 


0.5029882014 


0.4982256270 


43 


0.5076570564 


0.5029184361 


0.4982669896 


44 


0.5074804578 


0.5028518541 


0.4983064678 


45 


0.5073118214 


0.5027882423 


0.4983441875 


46 


0.5071506208 


0.5027274063 


0.4983802636 


47 


0.5069963746 


0.5026691682 


0.4984148013 


48 


0.5068486423 


0.5026133653 


0.4984478969 


49 


0.5067070199 


0.5025598478 


0.4984796388 



It is not hard to see that expected limit in each of the columns is 
1/2 as g — > oo. 

These two observations combined together give the factor 2 9+2n /n 29 
in the Weil-Petersson volume asymptotics. Similar to ([3]) it is natural 
to assume that the ratio j^j^m behaves like a n g bn , and it works. 
Moreover, it appears that a n = 2- 6 tt- 1 / 2 is independent of n and b n = 
n - 7/2. In the table below C 9yn = (2g)\ {^) 9 g n ~ 7 / 2 ^: 
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Vg,2/^g,2 


t/ „ in „ 


VgA/^gA 


A (~\ 

40 


1 m am o /ion 

1.019018429 


1.002495904 


0.9962430037 


0.9999695265 


A 1 

41 


1 m or /i^/ioo 

1.018547428 


1.002435270 


0.9963349432 


0.9999703519 


42 


1.018099193 


1.002377513 


0.9964224911 


r\ r\r\r\r\'~7A 1 O Af\ 

0.9999711349 


A O 

43 


1.017b 72110 


l.UUzo224ol 


0.9965059531 


0.9999718768 


A A 

44 


1.017264718 


1.002269843 


n r\r\n ro r^nno 

0.9965856093 


0.9999725812 


A C 

45 


l.Ulbo/oboo 


1.00221y5o4 


o.yyobbi ( i4o 


o.yyyy /325io 


46 


1.016503797 


1.002171501 


0.9967345009 


0.9999738890 


47 


1.016147948 


1.002125457 


0.9968041809 


0.9999744975 


48 


1.015807118 


1.002081325 


0.9968709494 


0.9999750781 


49 


1.015480379 


1.002038988 


0.9969349845 


0.9999756329 



We see that the ratio V^ n /C 9) „ apparently tends to 1 as g — > oo 
for any n = 1,2,3,4 (a standard extrapolation gives 1 up to at least 
6 decimal digits). For other values of n the situation is the same. It 
is worth mentioning that the case n = is computationally harder 
because of an additional non-trivial integration [7j, so currently we are 
able to compute V 9! q only up to g = 30. However, Conjecture 1 is 
rather precise even in this case: 
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21 


1.091195176 


22 


1.086790774 


23 


1.082792056 


24 


1.079145459 


25 


1.075806445 


26 


1.072737684 


27 


1.069907649 


28 


1.067289535 


29 


1.064860399 


30 


1.062600498 



Clearly, this sequence converges and its evaluated limit is 1 as well. 
Considerations that led us to Conjecture 2 are very similar to the 
ones described in this section. 
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